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Abstract 

We study the static phase structure of the ISS-KKLT model for moduli stabilisa- 
tion and uplifting to a zero cosmological constant. Since the supersymmetry breaking 
sector and the moduli sector are only gravitationally coupled, we expect negligible 
quantum effects of the modulus upon the ISS sector, and the other way around. 
Under this assumption, we show that the ISS fields end up in the metastable vacua. 
The reason is not only that it is thermally favoured (second order phase transition) 
compared to the phase transition towards the super symmetric vacua, but rather that 
the metastable vacua form before the supersymmetric ones. This nice feature is ex- 
clusively due to the presence of the KKLT sector. We also show that supergravity 
effects are negligible around the origin of the field space. Finally, we turn to the 
modulus sector and show that there is no destabilisation effect coming from the ISS 
sector. 
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1 Introduction 

In the last years, quite a large attention has been given to the problem of moduli 
stabilisation, especially concerning the cosmological implications they could have [lj. 
Following an earlier proposal [2j, Kachru, Kallosh, Linde and Trivedi (KKLT, [3]) 
have recently provided the first explicit model in which all moduli are fixed. They do 
so by turning on fluxes in a first step, which fix the complex moduli and the dilaton 
S, and introducing non-perturbative superpotentials [I] in a second step in order 
to stabilise the Kahler moduli T. For a more detailed study of the phenomenology 
arising from these models, see [5]. 

Unfortunately, the resulting low energy potential for T has an anti-de Sitter vac- 
uum which needs to be uplifted. The strategy proposed in [3] was to introduce an 
anti-L>3 brane far away from the visible sector in the sense of the compact dimensions 
so that the fine-tuning of the cosmological constant would be naturally explained by 
their curvature. However, such a mechanism results in a non-linearly realised super- 
symmetry, and therefore the low energy theory can no longer be expressed in terms of 
usual 4D supergravity. Latter attempts tried to realise the uplift by D-terms [6], i.e. 
using a fully supersymmetric sector, but this generically leads to a heavy gravitinc0- 

Parallel works have considered an F-term uplifting [;7j . This relies on adding a new 
sector in which supersymmetry is spontaneously broken by some field 3>, F$ ^ 0. If 

1 Actually this is not the case for the last two references of [6] because the uplift there is mainly realised 
by an i^-term. 
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this sector and the KKLT setup are decoupled in such a way that the Kahler potential 
and superpotentials add up, then the uplifting is trivially realised by relating the 
parameters of both sectors. The last two years, a rather large sample of such models 
together with their direct phenomenology have been proposed [8"j I§1 [TU]. 

In this paper, we shall focus on the setup developped in [SJ, where the uplifting 
sector was chosen to be the Intriligator, Seiberg and Shih model (ISS, [11]). A non- 
exhaustive list of string realizations of it can be found in Ref. [12]. This dual SQCD 
model is of particular interest since it realises a breaking of supersymmetry in local 
minima in the squarks direction. Elsewhere in the field space, in the mesons direction, 
there are super symmetric vacua, and both locally stable points are separated by 
a potential barrier. This ensures a long life-time for the SUSY breaking vacua. 
Hence, one does not have to give up the idea of a global (supersymmetric) minimum. 
Now, from cosmological considerations, we are led to wonder whether we ended 
up living in the metastable vacuum or not. And indeed, following the results of 
Refs. [1 3|. 1141 [151 I16j . we do. There, the authors showed that finite temperature 
correctional to the ISS setup favour the fields to go in the metastable vacua rather 
than in the supersymmetric ones. In [TJ], it was assumed that the fields start in 
the supersymmetric phase. Instead, the authors of [151 E] assumed the starting 
point to be the origin of the field space, which is a minimum at high temperaturqi 
This is because the origin of the field space contains the highest number of light 
degrees of freedom and hence maximises the entropy. We shall adopt the same 
attitude. In these various studies, it was found that the supersymmetric vacua form 
at a higher temperature than the metastable ones, but the origin is always a local 
minimum in the mesons direction. Therefore, the phase transition is first order 
towards the supersymmetric vacua. This is thermally disfavoured in comparison 
with the second order phase transition that occurs towards the non-supersymmetric 
vacua, even though the latter happens at a lower temperature. 

In this paper, we complete the study done in [16]. We work out the complete 
phase structure of the ISS-KKLT model. As will become clear in the text, the ISS 
fields do end up in the non-supersymmetric vacua, and the modulus, on the other 
hand, is not destabilised by thermal effects, as suggested in [18]. We will show that 
in our case, the presence of the modulus sector modifies the thermally corrected ISS 
picture in such a way that the metastable vacua form first, and they remain the true 
vacua of the theory during a certain time. Later on, the supersymmetric vacua form, 
but the fields have long gone in the SUSY breaking ones. Not until an even lower 
temperature are the two vacua degenerated. From that moment on, the fields can 
tunnel down from the metastable vacua to the supersymmetric vacua. 

Another feature that was pointed point out in [16] is that the origin of the ISS 
field space may no longer be a minimum at high temperature when this sector is 
coupled to the KKLT sector. This is due to supergravity, and could have a non- 
trivial effect on the phase transition. We study in great detail this point and find 
that as expected, this displacement is very small. 

However, let us emphasize that this study is still at the toy model level. We will 
not at all address cosmological problems such as the gravitino overproduction that 

2 For a review on finite temperature field theory, see [17j . 
3 We will develop these points in the following Sections. 
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usually happens when the super symmetry breaking sector is in thermal equilibrium. 
Even though the present paper obviously aims at a more realistic application, we 
leave these investigations for future work. 

The paper is organised as follows. Section [2] reviews the zero temperature ISS- 
KKLT setup in order for the paper to be self-contained. We introduce the main tools 
of finite temperature effective potential in Section In Section U the relevant tem- 
peratures and phase transition of the ISS sector are derived assuming the rigid limit 
(zeroth order in supergravity expansion). In subsection 14. 1( we compute the critical 
temperature of the second order phase transition towards the would-be metastable 
vacua. We give an insight of how the supersymmetric minima form in subsection 14.21 
We eventually compute the degeneracy temperature between the non-supersymmetric 
and the supersymmetric vacua in subsection 14.31 The rigid limit assumption of Sec- 
tion 2] is verified by working out the supergravity corrections around the origin in 
Section [5j Section [6] deals with the modulus sector. We show that the temperature 
corrections coming from the thermalised ISS sector do not destabilise the modulus. 
Finally, we conclude in Section [7] and draw the future directions that seem relevant 
to us. 



2 ISS-KKLT model 



Let us start by recalling the KKLT construction for moduli stabilisation in the frame- 
work of type IIB string theory. In [3], the authors used non trivial background fluxes, 
i.e. non-zero vacuum expectation values for certain field strengths in the internal di- 
rections, in order to stabilise all complex structure moduli as well as the dilaton. 
However, the Kahler modulus T, which describes the fluctuations of the overall in- 
ternal volume, cannot be stabilised in this manner. Non-perturbative effects such 
as gaugino condensation on D7 branes are used to generate an Affieck-Dine-Seiberg 
[3] superpotential at an intermediate scale A <C Mp. At low energy, the procedure 
results in the following setup 



K x = -3 In (T + T) 



Wi 



Wq + ae 



-bT 



(2.1) 



where the constant Wq is remnant of the stabilisation of all other moduli at the 
Planck scale. 

The model exhibits a supersymmetric minimum DtW\ = 8tW\ + K^W\ = at 
T = T , implying 



Wq 



a c 



-bT 



< , 



(2.2) 



KKLT/ 



K 1 



K tt Dt Wi D?pW i -3\Wi\ 2 



i 2 b 2 e- b ( To+T °) 
3 (T + T ) 



where K TT = [K 1 )- r ^ is the inverse metric for the Kahler potential K\. 



< , 



As mentionned in the Introduction, the energy can be uplifted to a positive value 
by adding a sector in which supersymmetry is spontaneously broken. In [8], the 
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uplifting sector was chosen to be the ISS model [TT] 

K 2 = Tr\cp\ 2 + Tr\p\ 2 + Tr|$| 2 , W 2 = h TV (<p$><p) — hfi 2 Tr $ . (2.3) 

This is the magnetic dual of a SUSY-QCD theory with gauge group SU(N C ). When 
the number of flavours satisfies Nf 3N c /2, the electric theory is asymptotically 
free whereas its dual, with gauge group SU(Nf — N c ), is infrared free. 

The magnetic fields under consideration are the gauge singlets <3? = f $j J , which 
we call mesons because they are in one-to-one correspondence with the electric 
mesons. The quarks cp = [fa), and the anti-quarks p = (iff) are in the funda- 
mental and antifundamental representations of SU(N). In the rigid supersymmetry 
limit, the theory has a global symmetry G = SU(N f ) L x SU(N f ) R x U(T) B x 
U(l)' x U(1) R which is explicitly broken to SU(N f ) x U(1) B X U(1) R by the mass 
parameter \i. 

We denote by N the magnetic number of colours N = Nf — N c , which satisfies 
Nf ^ 3iV. The indices run as i,j = 1, . . . , Nf and a = 1, . . . , N . For convenience, we 
will omit the flavour and colour indices from here on and will just keep in mind that 
$ is an Nf x Nf matrix, whereas if and ip^ are Nf x N matrices. 

The setup (|2.3[) has supersymmetry breaking solutions 

V = t T = (^) , * = (2.4) 

generated by non- vanishing F-terms for the mesons F$ = h ({pip — fj? Ijv ) • Notice 
that the supersymmetry breaking does not affect the gauge sector since it is driven 
by gauge singlets. The corresponding vacuum energy is 

V min = \h 2 ^\(N f -N) . (2.5) 

Far away from the origin in the mesons direction, after integrating out the quarks, 
gaugino condensation produces a non-perturbative superpotential [1] 

fh N f det$\ , \ 

which gives rise to super symmetric vacua 

(h$) = A m e 2N /( N f- N H Nf . (2.7) 

In the above expressions, A m is the dynamical scale of the magnetic theory, and 
e = ///A m is a small parameter. The existence of these vacua renders the non- 
super symmetric ones (|2.4p met ast able. Both regions of the ISS field space are sep- 
arated by a potential barrier. The lifetime of the metastable vacua can be made 
arbitrarily large by tuning e very small, or equivalently, A m very large for [i fixed. 

We now couple both sectors in the following way 

K = K X {T,T) + K 2 (x i ,Xj) , W = Wx{T) + W 2 ( X i ) , (2-8) 



4 



where \ % denote collectively the ISS fields ip, tp, <£. 

As explained in [8], such a decoupling between the two sectors can be achieved 
by considering systems of D3 and D7 branes. The gauge sector SU(N) arises from 
a stack of N D3 branes. Therefore the ISS gauge coupling, the dynamical scale A m 
and the mass parameter \x depend on the dilaton, which was already stabilised at 
higher energies. The mesons are interpreted as the positions of Nf D7 branes, this 
ensures the decoupling in the Kahler potential (|2.8p . The (anti-) quarks, on the other 
hand, are seen as open strings in the D3-D7 sector. Thus, their kinetic terms may 
not be canonical, but modifying the Kahler potential does not affect the main picture 
of the model since they are not directly related to the supersymmetry breaking, and 
hence to the uplifting mechanism. 

The supergravity corrections are negligible around the metastable vacua, as they 
are higher order terms in powers of the ISS fields ~ fi 2 /Mp. There, the scalar 
potential is well approximated by 

V( X \Xi,T,T) ~ 1 Viss{x\Xi) + Vkklt {T,T) , (2.9) 
(T + T) 

where VJss is the global supersymmetric (as opposed to supergravity) scalar poten- 
tial for the ISS sector. However, when computing the critical temperature, we will 
consider the expansion (|2.9p to be valid at the origin of the field space as well. We 
then explicitly verify it in Section [5j 

The fine-tuning of the cosmological constant to zero is given by 

(V) = =► |/iV|(iV/- N) ~ 3 1 Wol 2 , (2.10) 

and illustrated in Figure [TJ 

V(TW ^ n 




Figure 1: The KKLT potential (purple, dashed) and the uplifted ISS-KKLT potential 
(green, plain). The vev T ~ Tq is not modified by the uplifting mechanism. 



On the other hand, the gravitino mass is 



where (|2.2p together with the condition VTq ^> 1 were used in the last 



Here and in the following numerical results, we fix a = h = 1 and b = 0.3. We 
also set Nf = 7 and N = 2. Asking for a TeV range gravitino mass and imposing 
(|2.10p . the parameters are found to be 

T ~ 110 , |W | ^ (10~ 14 - 10~ 13 ) Mp , n ~ 1 - 5 • 10" 7 M P . (2.12) 

When needed, we will also take a coupling constant g = 0.1. Since the ISS gauge 
sector lives on D3 branes, one should have g{Mp) ~ 4-7r/ReS' and run this value 
down to a scale of order [i. However, the dilaton vev strongly depends on the UV 
completion of the model. Though we believe this point to be crucial, it goes far 
beyond the aim of the present work. We also consider the scale A m ~ 10 11 Mp. Let 
us recall that this scale is a Landau pole, which is not physical, and it is therefore not 
surprising to have A m higher than the Planck scale. The previous value corresponds 
to g (Mp) = 0.5 which means that the theory is perturbative at the Planck scale. 

Before coming to the specificities of our model, let us introduce the finite tem- 
perature formalism. 

3 Finite temperature effective potential 

The general one-loop effective potential including finite temperature effects can be 
split into different contributions [19j 

V eS (x\T) = V (x\T) + V?{ X \T) + V®{ X \T) , (3.1) 

where 

V = e K {k tt D t WD t W + K fj DiWDjW -3\W\ 2 } (3.2) 

is the tree-level supergravity potential, and K tJ is the inverse metric for the Kahler 
potential K2. 

The potential V® is the usual one-loop temperature independent Coleman- Weinberg 
effective potential [20] , and V® is the finite temperature contribution 

V e = ^_lj2nBj°°dx x 2 In (l - e -V* 2 +^/e^ 
KB ^ 

- E n F dx x2 ln i 1 + e-^ 2+M £ /e2 ) J . (3.3) 

Here rip (^_f) are the bosonic (fermionic) degrees of freedom, and Mp (Mp) are the 
bosonic (fermionic) field-dependent mass eigenvalues. One can immediately see that 
finite temperature corrections break supersymmetry. 

The potential (|3.3|) may be expanded at high temperature, 2 S> M|, Mp, 

^"^T ( nB + l nF ) + I" (3TrA^ + TrA4 + TrA^)+... , (3.4) 

4 As usual, we assume that the uplift of the modulus potential does not substantially modify its 
This can be easily verified graphically (Figure [lj. 
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where M x are the mass matrices for vectors, fermions, and scalars, expressed in 
terms of the fields. The trace Tr Ad'j is summed over Weyl fermions. 

In general, one should use the following supergravity formulae for the mass ma- 
trices in the presence of a non-canonical Kahler potential [2T] 



Tr M) = (e G \k ab K cd (V a G c + G A G C ) ( + -2\) , (3.5) 



and 

TrM 2 s = 2(K AB d ] Vo _ ) . (3.6) 
d x A dx B 

In the above expressions, x A represent the scalar fields in thermal equilibrium. The 
term —2 entering (|3.5p takes into account the mixed Goldstino-gravitino contribution 
|21j . The function G = K+hi \ W^ is the supergravity Kahler invariant potential, and 
we also introduced Ga = dG/dx A and V aGb = Gab ~^ab^c, with the connection 

T% = K CD "d A K B n . (3.7) 

However, as briefly mentionned in the Introduction, we will be concerned with 
these general results in Section [5] when we explicitly calculate how supergravity to- 
gether with temperature effects displace the minimum from the origin in the mesons 
direction, and in Section [6] when we study the destabilisation of the modulus. Since 
both T and <J> are singlets under SU(N), the gauge bosons contribution will not be 
relevant when computing the derivatives of the effective potential (|3.4p . This is why 
we did not write Tr M 2 here above. 

At the origin, we keep the ISS sector at the rigid level. When computing the finite 
temperature corrections there, we shall use the results of global supersymmetry 

3Tr M 2 V = 6(DfD ai ) , 
TiM} = (F^Fij) + 4(D?D ai ) , (3.8) 
TiM 2 s = 2{F ij F ij ) + 2{DfD ai ) , 

where, as usual, Fij = d 2 W/d\ l dx^ and Df = dD a /dx l - Here again, x % represent 
the scalar fields associated with ip, <p and <!>. One should be aware that the traces 
above run over the flavour and colour indices as well. The index a labels the adjoint 
representation of SU(N). 

Let us now turn to the main part of this paper, namely the phase structure of 
the model (|2.8I) once finite temperature corrections are included. 



4 Critical temperature and phase transitions 
in the ISS sector 

When the Universe cools down, the ISS fields end up in the non-supersymmetric 
vacua, as studied in [131 El EH USj - In this section, we show that the picture is 
not drastically modified when we add the modulus sector. However, this scheme 
is valid only if we consider the KKLT sector to be classical, which means that we 
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assume the modulus to be already lying in its minimum T = Tq and we neglect its 
quantum corrections to the ISS sector. In turn, Section [6] deals with the eventuality 
of a modulus destabilisation by temperature. 

Finite temperature effects are to restore all symmetries. At sufficiently high 
temperature, all the fields sit at the origin of the ISS field space. As we shall see, 
when the temperature lowers, the potential starts to exhibit a tachyonic direction 
towards the non-supersymmetric vacua, which form first. At a lower temperature, 
the would-be supersymmetric vacua form, but the origin remains a local minimum 
in the mesons direction (saddle point). Therefore, the origin and these new minima 
are separated by a barrier. 



4.1 Critical temperature 

In what follows, we focus on the behaviour of the potential at the origin of the field 
space. The symmetry restoration due to finite temperature appears when the tachy- 
onic tree level masses are compensated by the thermal masses (second derivatives 
of the potential (|3.3p ) at the origin. This is also a good reason to keep the lowest 
order in supergravity around the origin. Indeed, even though the corrections to the 
potential are negligible, supergravity effects could have a non-trivial impact on its 
derivatives and one typically has to take them into account. However, in the case of 
the critical temperature, and motivated by the results of [16], even if the exact loca- 
tion of the origin may vary with temperature and supergravity, the moment when the 
curvature of the potential at the origin becomes negative should not be drastically 
affected by supergravity effects. This, obviously, assumes that the origin is indeed a 
minimum at high temperature, as we will show in Section 

We use the high temperature expansion (|3.4p because the tree level masses are 
of order h 2 /j, 2 / (Tq + To) , see for instance fl2. 12f) . We follow the standard procedure 



|19j but there is no need to shift the fields here since we work at the origin of the 
field space. 

The traces (|3.8p expressed in terms of the fields are easily calculated from the 
superpotential (12, 3p . We find 

„ T KA 2 2 ^-1 TrM 2 + Tr|yf 

3TrA4„ = 3g — _ 3 , (4.1) 

N {To + To) 

. ,9 fr-? -?N 2 - 1\ TrM 2 + Tr|(o| 2 Tr|$| 2 
TvM 2 f = 2\h 2 N f +g 2 — -— — ^- + 2h 2 N U ^r , 



(To + To)" (To + Tq) 

where g is the coupling constant. It follows that the potential (|3.4p reads 



,e © 2 IV,, „ ,N 2 -l 



Vr= ~ _^ \(h Z N f + g''——\ TrM 2 + Tr|£| 2 
4 (To + Tq) 3 IV / 

(4.2) 

where we dropped the constant term cx 4 in (|3.4p since it is not relevant for the 
computation of the critical temperature. 



+ h 2 NTr\<S>\ 2 
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We compare the scalar thermal masses in (|4.2p to the tree level masses at the 
origin. The latter are ± h 2 /j, 2 / (Tq + To) or ± h 2 [i* 2 / (Tq + To) for the squarks and 
for the mesons, as easily seen from the superpotential (12. 3p and from the expansion 
(|2,9p . The thermal mass matrix is diagonal and positive definite, while the classical 
mass matrix is anti-diagonal. We ask for the determinant of the whole squared mass 
matrix to be zero at the critical temperature O c . This means that all the eigenvalues 
are positive above the critical temperature, while below C , tachyonic directions 
appear in the potential towards the would-be metastable vacua. From (|4.2p . we get 

which is in agreement with the critical temperature derived in [16J in the rigid limit. 
As computed there, this is only slightly modified by supergravity corrections. 

The critical temperature (|4.3p is of order fj 2 and the tree-level masses are of 
order hfs 2 / (Tq + To) , with Tq ~ 110. Therefore, the use of the high temperature 



expansion (|3.4p is a posteriori justified. 

For the values of the parameters given in (|2.12p . i.e. for a gravitino mass TO3/2 = 1 
TeV, we find that C ~ 4 • 10 -7 Mp. Notice that the critical temperature does not 
depend on the modulus, as one could have expected from the expansion (|2.9p . This 
is a crucial point since it is the reason why the would-be metastable vacua form first, 
as we show now. 

4.2 Formation of the supersymmetric vacua 

Having computed the critical temperature does not yet ensure that the ISS fields 
actually go in the metastable vacua. In this section, we turn to the mesons direction 
and work out the temperature SU sy at which the SUSY preserving vacua appear. In 
particular, we want to know if they are already formed when G = G c . In order to 
achieve this, we ask for the mesons to be away from the origin and integrate out the 
heavy quarks. The low energy theory is then pure Yang-Mills, it is strongly coupled 
in the IR and gaugino condensation [1] produces the non-perturbative term (|2.6p 
which gives rise to the supersymmetric vacua 

W NP = NA(det$) 1/N , (4.4) 

with A = h^A"^ 3 and v = Nf/N. 

At zero temperature, the vacua are $0 = (A~~ 1 h[i 2 ) 1 ^ u ^ Ijv, and the quarks 

3/2 

masses are = h$o/ {To + To) 

We decompose the mesons into a classical background $ and a quantum field 
as follows 

$ = ®l Nf + <j) . (4.5) 
Expanding the total superpotential /iTr (tpQip) — hfi 2 Ti + Wnp according to the 
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above decomposition, we get 



W = 



+ hTr (jjxfxp) + i A$ 



NA&- 1 - hfi 2 N f ) $ + ( A$ 




Tr^ 2 



hfi 2 ) Tr<p + M>Tr (pip) 




(4.6) 



Keeping the quadratic order in <p is sufficient because, following the standard proce- 
dure [19], we express the masses in terms of the classical field <3? and hence higher 
powers in <p are not relevant. Notice also that the quarks ip and <p should not be 
present in W since they have been integrated out. This point will become clear as 
we advance in the computation. 

At this stage, we would like to emphasize that working out the whole finite tem- 
perature corrected potential in the context of supergravity drives a lot of technical 
complications. For the sake of clarity, in order to sketch the mechanism that hap- 
pens around the SUSY vacua, we will again consider the rigid limit. We believe that 
supergravity corrections do not strongly modify the following results. 

If the quarks are integrated out, it means that the temperatures we consider are 

3/2 

6 <C h$>o/ {To + To) . On the other hand, using (|3.8p . we find that the masses of 
the mesons are 



which is in agreement with [15] . The high temperature expansion (|3.4p is thus legit- 
imate and one finds that the effective potential is 



From this expression, we see that the thermal contribution to the mesonic mass 
matrix, namely the second derivative of (|4.8p with respect to <1>, is diagonal and 
positive definite. Therefore there is no way that this contribution can lead to a 
destruction of the SUSY vacua and hence to a "critical" temperature. In other 
words, the masses are already positive at zero temperature and thus the origin and 
the supersymmetric vacua are separated by a barrier. When temperature effects are 
included, only a first order phase transition can happen. 

Moreover, the contribution (|4.8|) is much smaller than that of the tree-level masses 
rrifj, ~ A<& v ~ 2 1 (To + To) 3 ^ 2 by assumption^. Even if we used a high temperature ex- 
pansion, recall that the thermal mass is proportional to the quartic self-coupling of 
the mesons. Since it arises from non-perturbative effects, it is unnaturally small. 

Even though the squarks have been integrated out at tree level, their effect in the 
loops may be important. However, due to their large masses, the high temperature 

5 Since T ~ 110, the assumption 6 <C h$j (T + T ) 3 ^ 2 implies that 9 is even smaller than h$. 




(4.7) 



{To + T ) 
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expansion can not be used. From (|3.3p . we derive a low temperature expansion 



9 5 / 2 



3/2 -M B /6 



+ J2n F Ml /2 e 



3/2 -M F /Q 



(4.9) 



with M Bj f > 6. 

At the leading order, the squarks mass matrix is almost diagonal and its eigen- 

g 

values are h 2 & 2 / (Tq + To) . Recall this sector is super symmetric, so the fermionic 
degrees 

3NN f 



and bosonic degrees of freedom give the same contribution. Using (14, 9j) . we find 

\ 3 / 2 ( 



e 



7T' 



3/2 



{To + T ) 



3/2 



05/2 



cxp 



e (T + T ) 



3/2 



(4.10) 



This consists of a negative contribution which, together with the tree level potential 
deduced from (|4.6|) and with the effective potential (|4.8|) . gives the total potential. 
One then has to consider the system 

dv tot /d$ = o = dV tot /d$ 2 . 

Solving it brings us to knowing the temperature SU sy and the corresponding vev 
$ (©susy)- However, it turns out to be very hard to solve and we approximate = <I>o 
at all temperatures. We concentrate on the second equation of the expression above 
and find ^ 

tS 
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with 



B 



In 



N 2 



susy 



2v 



V 2 I 77' 



B (T 

3/2 



T f 2 



(4.11) 



A 2 $g"- 6 (To + T ) 3 



h^NNj 



> 



For B > 1, we fulfill the consistency condition that the squarks are integrated out 
at tree level. 



Numerical results 

As we already explained at the end of Section [21 the dynamical scale A m of the theory 
relies on the UV completion of our model. Since it is not the object of this work, we 
choose to consider the case of a half-unit gauge coupling at the Planck scale. Then 
A m = M P e- 2n ^ 3N - N f^ 2( - Mp ^ is approximately 10 n M P . Recall that N f > 3iV so 
that the argument of the exponential is positive. 

For this value and the rest of the parameters given by (|2,12p , we find the following 
results 

$o - 1-KT 3 , m Vt $ ~ 4-KT 7 , m ~ 3 • 1(T 13 (4.12) 
for the vev of the mesons and for the tree level masses, and 

B ~ 15 , 6 susy ~ 3-10- 8 
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for the temperature. All these results, except for B, are expressed in units of Mp. 
Notice that B is larger than one. 

The main conclusion is that the SUSY vacua form at a temperature which is 
smaller than the critical temperature (|4.3|) . Obviously, this result depends on the 
choice of A m and we emphasize, again, that a closer study of the UV physics of our 
model is required. However, we find numerically that the constant B is negative for 
A m smaller than 10 4 which corresponds to a gauge coupling of 0.9 at the Planck 
scale. Therefore it seems that the more the theory is perturbative at high energy, 
the more consistent the picture is. Varying A m in this range yields 

10 4 < A m < 10 15 => 1.5 -10~ 8 < m<p,£ < 7.4 • 10" 6 . (4.13) 

Hence, even for very high values of A m , the upper bound under which one can 
integrate out the squarks is only slightly above the critical temperature (|4.3f) . and 
the corresponding SUSY temperature is 3 • 10~ 7 Mp < C . 

It is clear that the major cause of such an effect is the explicit dependence of the 
SUSY temperature on the modulus. This pushes the tree level squarks masses to 
very low values compared to the original ISS scenario. 

Another result that we were able to derive numerically is that already once 
the squarks are integrated out, their contribution (14.10p is very small compared 
to the tree- level one (|4.6p . This means that whenever one can consider the non- 
perturbatively generated superpotential (14. 4h . then the vacua are already there. As 
such, the SUSY temperature (14. lip does not really make sense, and we are more 
encline to rely on the evaluation of the squarks masses fn^s as i n (|4.13p . Also, since 
these are tree level masses, they do not depend on <3? (0) and thus are not biased by 
our approximations. 

The conclusion is unchanged : the supersymmetric vacua form after the would-be 
metastable ones, and this is due to the presence of the modulus. 



4.3 Degeneracy between the vacua 



Finally, in this paragraph we compute the degeneracy temperature ©deg a t which it 
becomes possible for the fields to go from the metastable vacua to the supersymmetric 
ones. This temperature is defined as the moment when both vacua have the same 
energy. 

The total number of degrees of freedom in the non-supersymmetric vacua is 
(Nf + N) 2 — 1. The vacuum energy there is 



(V)\ 



mcta 



vr 2 e 4 

24 



(N f + Nf 



+ 



/iV| {N f - N) 
(Tq + Tq) 3 



where we did not account for the KKLT energy since it is constant over the whole 
ISS field space. 

In the last paragraph, we showed that the squarks can be totally neglected in 
the supersymmetric vacua. Therefore only the finite (high) temperature correction 
coming from the mesons is relevant. Recalling that these vacua have zero energy at 
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tree- level, using (|4.8p and assuming again that <3? = <3?o, one finds 

_ vr 2 6 4 - ~ ^ + A 2 $ 2 - 4 9 2 

^)U--^/+ 8 (T 0+ T ) 3 • 

Actually, it is easily seen from our numerical results (|4.12j) that the last term in the 
above expression is negligible. To good approximation, the degeneracy temperature 
is thus given by 



deg ~ y vr 2 (2NN f + N 2 -1) (T + f)3/2 ' ( - 1 ^ 

Using our parameters, one finds numerically Odeg — 7 • 10 -9 Mp. As before, the 
degeneracy temperature explicitly depends on the modulus, reason why it is so low 
compared to the critical temperature. We believe that this is a major improvement 
over the case of an isolated ISS sector. As was already noted in [Hj, the presence of the 
modulus enhances the lifetime of the ISS metastable vacua. We confirm this result 
here by showing that the supersymmetric vacua actually become the true vacua of 
the theory only at relatively late times. 



5 Supergravity and finite temperature correc- 
tions at the origin 

The computations of Section [5] have assumed that the origin of the ISS field space is a 
minimum of the potential at high temperature. However, as pointed out in [16], this 
is not as straightforward once one includes supergravity. Consider for instance the 
cross term K TT K T W 2 d T W l in the supergravity potential (|3.2j) . It contains a linear 
term in $ which contributes as a constant to the equation d$V = 0, and produces 
a displacement from the origin. Generical temperature corrections contain similar 
terms and one has to work out the full supergravity plus temperature corrected 
potential and solve for a minimum around the origin. This is an important point 
because, even though unexpected, the displacement could be large enough to spoil 
the phase transition towards the supersymmetry breaking vacua. 

From the superpotential and Kahler potential ([2.3)1 ; one can see that only terms 
of at least quadratic order ~ y? 2 , <^ 2 , iptp can appear in the scalar potential. Con- 
sequently, the origin (p = (p = is always a solution to the extremum equations 
d^V = = d^V. In what follows, we concentrate on the equation d$V c g = in the 
background ip = (p = (here V e g stands for the full potential defined in (|3.ip ). 

6 This result is in slight disagreement with [To]. First of all, they computed the degeneracy temperature 
from the origin to the supersymmetric vacua. Indeed, when the KKLT sector is not present, the latter 
form before the metastable vacua. However, by dropping the T-dependence and replacing Nf by Nf — N 
in the prefactor of (|4. 14|) . we do not find exactly their result. This is due to the fact that they use a high 
temperature expansion even for the squarks in the supersymmetric vacua, which results in dropping the 
2NNf in the denominator. 
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The tree-level scalar potential 



Vn 



K TT D T WD T W + K'WiWDiW - 3 |Wf 



(5.1) 



receives temperature corrections given by (|3.4p . where the mass matrices squared 
(j3.5p and (|3.6p can be developped using the semi-canonical Kahler potential f|2.8j) 

TiM 2 f = { e G \K^Ki\Gij + dGj) {G M +G- k G T ) - 2 

d 2 V , 



and 



(2RV- 



dx i dx j 

The new minimum at high temperature satisfies 
dV n 9 2 d 



d$ + 24 d$ 



{TrMj + TrMt} 



(5.2) 
(5.3) 

(5.4) 



93=1,0=0 



Let us start with the zero-temperature potential (|5.1|) . Differentiating with re- 
spect to <3? yields 



K TT {(D T WK<$, + Wq,K T ) W T + K^W T K T W) 
+ K fj DjW (K$DiW + KiWs, + W w ) + WD^W 



Since we expect the displacement (3>) to be small, it is sufficient to keep the linear 
order in $. One gets 



(d*V ) 



$ \k tt {DtWiW^ + WtK^Wi) + |/»V| JV> + l^il 2 } 
- hfi 2 l Nf {k tT K t W t + Wi- 2hpL* 2 Tr ¥} 



(5.5) 



where K\ and W\ are the pure KKLT potentials defined in (|2.1|) . 

It is a long but straightforward computation to derive the other two contributions 
in (|5.4p ; some steps are given in the Appendix [A] for the interested reader. 

The general solution to the linearised equation (|5.4p is of the form 



$ (T, T, 6) 



, 2I A + BQ 2 \ 



(5.6) 



where A,B,C,D are functions of T and T only, and given in (IA.6p . 

Figure [2] shows the behaviour of (|5,6p with respect to temperature for T = Tq. It 
is of some relevance to consider two different situations. For instance, the gravitino 
mass (|2.1ip fixes all the parameters, since the relation (|2.2p between To and Wq 
on the one hand, and the zero cosmological constant (|2.10p on the other hand are 
conditions of our model. 

We choose to consider m 3 / 2 = 1 TeV (blue, dashed line) and m 3 / 2 = 100 GeV (red, 
plain line) as an example. In both cases, as expected, the origin is the only vacuum 
at very high temperature. One can already approximate <I> ~ at © ~ 10~ 6 Mp for 
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3>(0) 

3>(0 C ) = 2.4 10" 1: 
3>(0 C ) = 0.2 10" 1: 




C =1.1O" 7 C =4.1O" 7 



Figure 2: Evolution of the mesons vev with the temperature for two values of the gravitino 
mass. When the temperature hits its critical value, the minimum at $ turns into a saddle 
point. 



the light gravitino case. The surprise comes from the fact that the minimum fades 
away from the origin very fast when the temperature lowers down, and this happens 
while the high temperature expansion is still valid. However it could very well be 
that keeping the linear order in $ is no longer a good approximation there. 

The critical temperature (|4,3p developped in Paragraph 14.11 depends on fx and 
thus on the gravitino mass. From Figure it is clear that the mesons still have a 
very small value at the critical temperature, for both cases we considered. Therefore, 
the phase transition towards the non-super symmetric vacua, in the squarks direction, 
will not be affected by the displacement. A problem would have arised if the mesons 
vev had been too high (at the critical temperature), forcing us to take into account 
the non-perturbative superpotential (|4.4p . 

Moreover, the system with a light gravitino remains around the origin during 
a longer time, ensuring even more the phase transition. Indeed, one could find a 
set of parameters matching our two conditions (existence of a minimum for T, zero 
cosmological constant) for any gravitino mass. In the case of a substantially heavier 
gravitino, not only the volume modulus would have a too small vev, but the phase 
transition towards the would-be metastable vacua would be spoiled. We conclude 
that, even though it is not a very strong effect, our model seems to prefer a light 
gravitino. 

6 The modulus sector 

Up to now, we have been considering that the modulus T was sitting in its minimum. 
In this section, we shall derive the condition under which this is valid at the typical 
ISS temperatures. 

Let us recall that since the moduli are only gravitationally coupled to the thermal 
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bath, their interaction rate is 



g 3 TT e 2 



As such, the moduli potential is not in thermal equilibrium. However, indirect tem- 
perature corrections coming from other sectors could destabilise a modulus because 
they would result in an extra source of uplifting. 

For instance, in [18] . the authors studied the maximal (or critical) temperature 
beyond which a minimum generated by non-perturbative effects would be destroyed. 
Assuming that the visible sector lives on D3-branes, the gauge coupling is directly 
related to the vev of the dilaton g 1 ~ 1/Re S. This implies that the dilaton potential 
is thermally perturbed through the gauge coupling. Typically, these effects destroy 
the minimum if they compensate the barrier between the metastable vacuum at 
Re S ~ 2 and the minimum at infinity. Therefore, in [18], the dilaton was destabilised 
for temperatures > ^Jm^Mp — 10 _8 Mp for a gravitino in the TeV range. The 
same could happen to the T-modulus if the visible sector lives on D7-branes. 

In a KKLT setup, however, the dilaton is not stabilised by non-perturbative effects 
W ~ e~ bS , but rather by non-trivial background fluxes W ~ m + nS. Whereas 
gaugino condensation takes place at a scale A <C Mp, resulting in a low mass for the 
modulus (T or S according to the model), a stabilisation by fluxes happens at high 
energy ~ Mp. The dilaton is then heavy enough not to be affected by temperature, 
and we can simply decouple it at low energy, as in the zero temperature theory. In 
what follows, we assume that this is the case, i.e. that the visible sector does live on 
D3-branes. 

In our model, the T modulus potential receives temperature corrections from 
the ISS sector. If there exists a critical temperature above which the potential is 
destabilised, we assume it to be higher than the temperatures computed in Section 
2J In this case, the ISS fields are at the origin, with the mesons slightly displaced, 
eq. 

We define the destabilisation temperature and the corresponding value T d for 
the modulus as the point where the minimum turns into a saddle point : 

^(T d ,O d ) = = ^f(T d ,Q d ) , (6.1) 

where V e s is the effective potential (|3.ip . 

The computation follows similar steps as in Section [5] and Appendix [Aj We 
simply give here the result for the effective potential at linear order in the mesons 
displacement 

Vq = Vkklt + e Kl [\h 2 ii*\N f -hiJ, 2 ATr$-hiJ.* 2 A* r Ir'§] , 

where A (T, T) = K TT K T W T + Wi was also defined in the Appendix El 
The traces of the mass matrices ( 15.21) and ( 15.31) are 

TiM} = -2e Kl 



|Wi | 2 - V 2 ^iTr$ - V* 2 ^iTr$ 



7 Since we expect the destabilisation temperature to be very high, the linear approximation made in 
Section [5] is even more valid as one can convince oneself from Figure [21 
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and 



TrMi 



2e 



A'l 



2N f (N f + 2N) {k^DtWiDtWi — 2 | Vl^i | 2 } 
+ {2 + 2N f (N f + 2N)} \h 2 fi 4 \N f 

'h^Ti $ {2N f {N f + 2N) + 1} (k tT K t W t + 2W X ) + h.c. 



This expression is easily implemented in a Mathematica routine in order to solve 
the system f)6. 1 1) . 

Veff(T) 



m 3/2 = 1 TeV 



100 110 120 T 

Figure 3: Non-destabilisation of the T modulus at high temperature. 



As a result, the one-loop effective potential is shown in Figured We point out that 
the constant term — C0 4 in (|3.4p has not been included for graphical convenience. 
One can see that there is no destabilisation of the modulus at all, and indeed the 
system (|6.1|) turns out to be non-solvable. 

We already argued in Section [5] that the parameters are fixed by the gravitino 
mass. In Figure El we took m 3/ / 2 = 1 TeV ; the rest of the parameters is the same, 
namely a = 1, h = 1, and b = 0.3. The constants Wq and jj, are fixed by equations 
(|2.2p and (|2.10p . and the solution $ (T, 0) was derived in the Appendix [AJ This 
computation assumes that the value of T at the minimum does not vary too much, 
which is cross-checked on Figure [3l 

In [22], the authors worked out the phase structure of the O'KKLT model [10J, 
which can be viewed as a simplified version of our model. Although they assumed the 
modulus to be in thermal equilibrium, it was found there that it is not destabilised 
by thermal corrections. In this perspective, we recover their result as the limit in 
which the thermal contribution of T is negligible, which is indeed the case of interest 
for an expanding Universe. 

7 Conclusions and future challenges 

Following earlier work [131 021 CE] i we have studied in great detail the static phase 
structure of the ISS-KKLT model when thermal corrections are considered. 
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We are now able to give the complete picture of its thermal evolution. At very 
high temperature, S> O c , the ISS fields are at the origin because this is the point 
where the entropy is maximised. At these temperatures, the modulus T is already 
stabilised (Fig. [3]). Once it lies in its minimum, we can consider it to be static and 
neglect its quantum corrections to the ISS fields. Then, as the Universe cools down, 
the ISS fields start being driven away from the origin (Fig. [2D, but they are still very 
close to it when the temperature hits its critical value C ~ 10~ 7 Mp. A second order 
phase transition takes place towards the would-be metastable vacua which at this 
stage are the true vacua of the theory. At a lower temperature, the supersymmetry 
preserving vacua form. They are separated from the origin by a barrier. Therefore, 
even if one enhances the dynamical scale A m in such a way that these vacua form 
first, it would consist of a first order phase transition and would thus be thermally 
disfavoured. For our parameters, however, the non-supersymmetric vacua form first. 
At a temperature 0deg ~ 10~ 8 Mp, the supersymmetric vacua become the global 
vacua of the model and from that moment on, the ISS fields can tunnel from the 
metastable vacua to the supersymmetric ones. 

Even though we tried to give a complete and quantitative study of the model, 
there are still challenges that deserve further attention. First of all, we showed that, 
if the visible sector lives on D3 branes, the modulus T is not destabilised by finite 
temperature corrections coming from the ISS sector. This assumes that the sector 
responsible for the stabilisation of T is out of thermal equilibrium. Another limiting 
point that we have not treated is the dynamical evolution of the system, especially in 
the modulus sector. Indeed, the potential generated for a modulus is generically so 
steep that it seems very unlikely that the field will actually end up in the minimum, 
and not overshoot the barrier towards the runaway minimuir|§ (this effect is known as 
the Brustein-Steinhardt problem [231 H]). Both issues have been recently addressed 
in (23J. Based on the conclusions of [IB], the authors have studied the conditions 
under which a stabilising sector (in their SUSY-QCD) in thermal equilibrium 

can lead to a destabilisation of the modulus at some temperature. They developped 
the whole set of dynamical equations when the stabilising sector is included in the 
thermal fluid, and constrained the initial conditions for the rolling modulus to reach 
its minimum. Their conclusion is that there is a region of initial conditions which lead 
to a stabilisation of the modulus. The allowed region is slightly reduced compared to 
the case where temperature corrections are not considered, but this is not a dramatic 
effect. We believe that these conclusions can be applied to our case - actually the 
authors of [23] do study the KKLT setup - knowing that on the other hand we 
have showed that the temperature contribution coming from the ISS sector does not 
destabilise T. However, we think that a closer evaluation of the dynamics of our 
model needs to be done. In particular, thermal fluctuations around the origin might 
be very important. 

Another interesting direction is inflation. It has been a big challenge for quite a 
while to combine inflation with string-inspired supergravity models : see for example 
[10], [25] and [26] for a review. Here, the coupling of the ISS-flaton [27] to supergravity 
as in the ISS-KKLT setup could be of particular interest [28J. 

8 I thank Z. Lalak and S. Pokorski for bringing my attention on this problem. 
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A Expression of the displacement of the mesons 

In this appendix, we derive the displacement <3? of the mesons in terms of the modulus 
T and the temperature G as given in (|5.6p • As already sketched in Section [5l we 
have to solve the equation 



^ + — 7^- TrM 2 f + TrVW< 







(A.l) 



<p=tp=0 



where Vq, Tr M. 2 and Tr.A/f 2 were respectively defined in (|5.ip . (|5.2p and (|5.3p . 

Keeping the linear order in it is easy to show that the tree-level (and thus 
temperature independent) contribution to the displacement is 



$ {k tt {D t W{W t +W t K t W)) + |/iV 4 | Nf + \W X \ 2 



hn 2 l Nf {, 



K^KtWt + Wi- 2 V* z Tr $ 



* 2itl, X. 



We now turn to the fermion mass matrix and compute <9$ 
term in (|5,2p gives the following contribution at the linear order 



TiM} 



( 3* { e G K fk K^ (Gy + GiGj) (G~ kI + G- k G T ) } 



V 

(2e Kl [¥{/i 2 (N + \fi 4 \N f )} + l Nf |^V|Tr$ 



The last term is 2e G which simply gives 



(A.2) 
The first 



(A.3) 



(5$ (-2e G )) = -(2e Kl $ \Wi\ 2 - hfi?l Nf {W x - V* 2 Tr *) ) 



(A.4) 



TiMj 



in equation (jA.ip . 



All together, (jA.3[) and (lA.4j) give the contribution <9$ 

The trace of the scalar mass matrix squared is given in (|5.3p and needs the same 
treatment as before : 



d 2 V _ 

dx l dx j 
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However, with some patience, one can get the following result for this contribution 



where we used the fact that K l] = 2Nf (Nf + 2N), which is a trace over the ISS 
scalar fields. 

From all these results, it is clear that the linearised solution takes the form <3? = 
$oIjVf which implies that l^Tr^ = A/$ . 

Eventually, plugging the different contributions into (jA.ip . the displacement of 
the mesons takes the form 



<2e 



¥| (4 + 2N f (N f + 2N)) ■ |/iV| N f + 2h 2 N 

+ (1 + 2JV>(JV> + 2N)) (rTTDtW^Wx - |Wi| 2 ) } 
hn 2 l Nf { - 6 (1 + N f (N f + 2 AT)) hfi* 2 Tr ¥ 

+ (1 + 2N f (N f + 2N)) (k tT K t D t W 1 - Wi) }] ) , 



(A.5) 




with the following entries 




(A.6) 



B (T, T) 
C(T,T) 



— {l + 2N f (N f + 2N)}(K TT K T D T W 1 -W^j 
K^lDrWtl 2 -2\W 1 \ 2 + 3\h 2 fi 4 \N f , 



Wi 



d It, t 



— [{1 + 2N f (N f + 2N)} {k tt \DrWil 2 - |W X 



) 



Wi 



2 



+ 3h 2 N + | /i V \N f (11 + 8N f [N f + 2N))] 
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